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Abstract. This work devotes to investigations of solar corona and solar wind plasma by the method of radio probing with using 
of very long baseline interferometry (VLBI). We performed the theoretical calculation of power spectrum of interferometric 
response to radio source emission, passed through the turbulent medium. Data of theoretical analysis are compared with results 
of international VLBI experiments on investigations of solar wind plasma. The observations were realized in 1998, 1999, 2000 
with participation of radio telescopes, included at Low Frequency VLBI Network (LFVN): Bear Lakes (RT-64, Russia), Puschino 
(RT-22, Russia), Urumqi (RT-25, China), Noto (RT-32, Italy), Shanghai (RT-25, China) and others. Preprocessing was carried 
out with using of S2 correlator at Penticton (Canada). Post processing of experimental data was performed at RRI (Russia) and 
was aimed to obtain value of solar wind velocity and index of spatial spectrum of electron density fluctuations. 



1. Introduction 

One of the most important problem of solar-terrestrial physics 
is the investigation of coronal plasma — the determination of 
physical parameters of electron density iiTegularities and their 
spatial-temporal structure. The method of radio sounding is 
widely applied for investigation of coronal plasma and solar 
wind. In this method the signal from the radio source goes 
through the turbulent medium and then is received at ground- 
based radiotelescope. Heterogeneous mediums cause phase, 
amplitude and frequency fluctuations of the received emission 
and distort the output signal of the instrument. The analysis of 
these disturbances allows one to get information on physical 
characteristics of the propagation medium. 

Perspective technique in this field proves to be inter- 
ferometry method [Altunin et all, 2000, Girin et all, 1999, 
Spangler S. R. et all, 2002], when sounding signals are prop- 
agated through circumsolar plasma by different paths and re- 
ceived at separated antennas of VLBI-complexes. A very long 
baseline interferometer permits studying the electron-density 
irregularities using relative phase, amplitude, and frequency 
fluctuations originating on propagation paths from the source 
to the receiving antennas of the interferometer Moreover the 
size of the baseline projection of the instrument to the wave 
front determines the scale of iiTegularities (100- 10000 km), in- 
terferometer being sensitive to them. Using of interferometric 
complexes with various orientation of baseline projections al- 
lows to obtain information about spatial structure of solar wind 
plasma parameters. The observations of sources, located at dif- 
ferent positional angles and distances from the Sun, let us make 
certain conclusions about the spatial structure and anisotropy 
of electron density of investigated medium. The radiointerfer- 



ometric method allows for obtaining data by receiving both 
monochromatic signals from spacecraft and wideband emis- 
sion from natural radio sources. 

In recent years, a number of VLBI-experiments connected 
with the investigation of solar wind plasma by the method of 
radio sounding has been carried out at Low Frequency VLBI 
Network (LFVN). 

We performed theoretical analysis of power spectrum of 
interferometric response to radio source emission, passed 
through the turbulent medium [Gavrilenko et all, 2002, Alimov 
et all, 2004] . The conclusions were compared with experimen- 
tal data having for an object to get information about physical 
characteristics of propagation medium, namely, solar wind ve- 
locity and spectral index. 

2. Frequency spectrum of the interferometer 
response to wideband emission under 
conditions of an inhomogeneous medium 

The purpose of current work is to obtain theoretical formu- 
las describing the response of an interferometer to the radio 
emission from the source, which propagated through a turbu- 
lent medium, to simulate VLBI experiments on radio raying of 
the solar wind, and to compare the calculation results with the 
experimental data. 

We analyzed the situation, when the noise emission of a 
cosmic radio source passes along axis z through the turbulent 
medium of solar wind with chaotic large scale irregularities of 
electron density. Then it is received by a two-element ground 
radio interferometer on plane (x, y), which is perpendicular to 
the direction of propagation. 
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The signals received at separated antennas having been 
transformed in the tracts of the interferometer, undergo mu- 
tual correlation processing. The power spectrum of signals at 
the output of the interferometer's correlator, having informa- 
tion about the propagation medium, can be represented in fol- 
lowing way: 



Y(Q.o 
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R(t) = {Ei(r, t) Elir + p, t) Ei(r, t + t) E*^{r + p,t + t)) (2) 

— the correlation function of the interferometer's output sig- 
nal, E\2 ~ e*' - — the signals at the outputs of the receiv- 
ing tracts of the interferometer, (D — the complex function 
of amplitude and phase fluctuations, caused by the turbulence, 
p — the baseline of the interferometer, Qo — the frequency of 
Fourier-analysis. 

Expression Q is described the correlation procedure, com- 
mon for interferometer reception. It has some advantages in 
comparison with traditional one point receiving as it allows to 
investigate the field fluctuations caused by turbulent medium 
only on two different propagation paths. In this case the influ- 
ence of source self-radiation is excluded and it allows us to 
sound the medium not only by monochromatic but also wide- 
band source signal. It proves to be very important when the 
later is raying by sound emission of natural sources. 

In this work the theoretical calculation of the power 
spectrum of the interferometer's output signal is performed. 
The distribution of field fluctuations was carried out with 
geometric-optical approach. It was supposed, that the tempo- 
ral changes of the electromagnetic field parameters in the solar 
corona concern the "frozen-in" hypothesis — when the irreg- 
ularities move from the Sun in radial direction at the velocity 
of solar wind V. The spatial spectrum of fluctuations of turbu- 
lence parameters F(k) was described by a power law function 
at the wave number range {k(),k„,], kq - In /A, K,n - 2nll,„ (A 
and Im — are outer and inner scale of turbulence respectively): 



F{k) = 0.033 cI{kI + 1?) 
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(3) 



^ — the structural coefficient, determining the intensity of 
fluctuations. 

The spectral index is usually considered to be equal to p - 
11/3 (Kolmogorov's spectrum) at distances from the Sun more 
than \5 ■ Rq (Rq — radius of the Sun) and it is decreased to 
p - 3 at shorter distances. Index p is the characteristics of the 
turbulent medium being investigated, which is determined from 
the spectrum of the interferometer's output signal. 

The expression for correlation function, considering the 
amplitude and phase fluctuations, has a form: 



oo 

■I 



R(t) = exp<^-4;rA^Z F{kJ[1 - cos k^V^t]x 



X [1 - cos{Kj_pj_)]dKj_ + i2n x 

Wo 



F{Ki_) smiK^Vi^r) sin(KiPj_) £/*:_l 



(4) 



A , <jj,, — plasma frequency, c — velocity of light, 

2Nuqc 

Ki., Pi., y 1. — projections of wavenumber, baseline and solar 
wind velociy on the plane of the sky, — the central fre- 
quency of the reception, — electron density, Z — the thick- 
ness of irregularities' layer. 

We examine two limit cases: strong and weak phase fluctu- 
ations. In case of strong phase fluctuations, the expression for 
y(Q()) can be written as follows: 
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In (|3-(|EJ the cross baseline projection is directed along 
the axis x; Jo, Ji — Bessel function. It is evident, that in case 
of strong phase fluctuations, the power spectrum of the output 
signal F(Qo) has the form of Gaussian function independently 
of the choice of the spatial spectrum's type F(k). Value (w^) 
described the intensity of phase fluctuations. 

The analysis confirms, that solar wind observations proved 
to be more informative at weak phase fluctuations. In this case 
the form the power spectrum must be drop-down power law 
function. When stream velocity is directed along baseline (p,,, 
Vy - 0) dependence of spectrum from baseline and solar 
wind velocity is described approximately by the expression: 
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(K — constant coefficient, (|9jl — is received under the condi- 
tion: Qo « InVxKm)- 

The function describing the spectrum is represented in fig. 1 
in logarithmic scale; the dimensionless value Q = Qo/Qi. is 
marked along the abscissa axis (Q,, - InVxIpx)- 

It follows from (|5Jl that the spectrum should have oscilla- 

'QoPx 



tions determined by the factor 



1 - cos 



and dependent 



only on the relationship between the solar-wind velocity and 
the interferometerbaseline. This makes it possible to determine 
the velocity of irregularity transport on the sounding path. 

Two characteristic intervals can be distinguished for this 
curve. The corresponding frequency dependences are described 
by the next approximation: 



< Qo < nVxIpx 



(10) 
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Fig. 1. The power spectrum at weak phase fluctuations 
F(Q)«[Qr^+i, Qo>7ry,/p,. 



(11) 



Therefore, the index p of the spatial spectrum of fluctua- 
tions of medium parameters can be determined from the slopes 
of the curve. Approximating curves ilOi and (lll> intersect at 
some point £1, from which one can determine the solar-wind 
velocity if p , is known. 

If the interferometer baseline is perpendicular to the drift 
velocity of irregularities of the medium (p,, » Pi), then 
yields 



F(Q) X const, < Q() < ttV^/Px, 
y(Q)«[QF+i, f2o>7ry,/p,. 



(12) 
(13) 



In this case, as well, the spectral index p can also be de- 
termined from the asymptotic dependence of the measured fre- 
quency spectrum, while the solar wind velocity can be deter- 
mined from the characteristic frequency at the break point of 
measured spectrum - jrVJlp^. 

3. Results of experiments 

We carried out the analysis of international VLBI experi- 
ments on investigations of solar wind plasma, implemented 
in 1998, 1999, 2000. These observations were realized with 
participation of radio telescopes, included at Low Frequency 
VLBI Network: Bear Lakes (RT-64, Russia), Puschino (RT- 
22, Russia), Hartebeesthock (RT-25, South Africa), Arecibo 
(RT-305, USA), Urumqi (RT-25, China), Noto (RT-32, Italy), 
Shanghai (RT-25, China), GMRT (RT-45, India). At the course 
of these experiments VLBI complex received the extragalactic 
source radio emission (wavelength 18 cm) passed through the 
turbulent plasma of solar wind. Radio sources were located at 
different angular distances from the Sun (4-37 degrees). The 
observational information were recorded using the S2 system. 
Preprocessing of VLBI-data was carried out at S2-correlator 
at Penticton (Canada). The data were calibrated in the NRAO 
AIPS package using the standard technique. 
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Fig. 2. The examples of the autocorrelation function (a) 
and power spectrum (b) taken in the course of experiments 
INTAS99.4 (November, 1999, source NRAO530, solar elonga- 
tion 6 = 16°, baseline Bear Lakes (RT-64, Russia) — Noto 
(RT-32, Italy), A - cm, integration time 0.1 sec.) 



Postprocessing was performed at RRI (Russia) and was 
aimed to obtain value of solar wind velocity and index of spa- 
tial spectrum of electron density fluctuations. The autocorrela- 
tion function of output signal was counted. After that Fourie- 
analysis was performed; as a result we have the power spectrum 
of field of interferometer signal. The examples of the autocorre- 
lation function and spectrum taken in the course of experiments 
are shown at fig. 2. 

The results of spectral analysis of realization series were 
compared with theoretical conclusions. Obtained experimental 
power spectra let us to get estimations of index of spatial spec- 
trum p according slope of the spectrums. During observations 
p took on a value 3.2-4.5 depending on the angular distances 
of radio sources to the Sun. 

The determination of solar wind velocity were embarrassed 
by the several reasons. The most part of experimental data was 
correlated with too large integration time (f = 2 s). This circum- 
stance did not allow to get output spectrum in sufficient wide 
band, where desired feature of spectrum is located. Reiteration 
of preprocessing with integration time f = 0.1 s, before activ- 
ity correlator was stopped, has been implemented for sources, 
observed at experiment INTAS99.4 in 1999, November. 

Strong response was obtained only for one source 
NRAO530, located at angular distance from the Sun 6 - 16 de- 
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Fig. 3. The power spectrum, obtained for source NRAO530 
(experiment INTAS99.4, 1999, = 16°, baseline Bear Lakes 
(RT-64, Russia) — Nolo (RT-32, Italy), ^ = 18 cm, integration 
time 0.1 sec); spectral index, evaluated from spectrum slop p - 
3.8. 



Experimental works are satisfactory corresponded with 
conclusions of theoretical analysis as regards to determination 
of spectral index. Nevertheless determination of average ve- 
locity of irregularities transportation is embarrassed by weak 
source signal and existing sufficient fluctuations of velocity. 
With this object the program of solar wind observations was in- 
cluded at schedules of last experiments VLBR04. 1 , VLBR04.2, 
VLBR04.3, performed at Low Frequency VLBI Network in 
June, July and October of 2004 at wavelength 6 cm. Records of 
experimental data was realized at system MK2. Processing of 
date will be made at correlator MK2 "NIRFI-3" at RRI, Nizhnij 
Novgorod, Russia. The high time resolution of this correlator 
and large diameter of enabled antennas let us to expect satisfac- 
tory results as concerned the solar wind velocity determination. 
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grees. Positions of other sources, coming for solar wind obser- 
vations, were too close to Sun and their signal-to-noise relation 
was too low, that does not allow to obtain distinct spectrums. 
The sample of power spectrum for source NRAO530 is dis- 
played on fig. 3 (baseline Bear Lakes (RT-64, Russia) — Nolo 
(RT-32, Italy), wavelength 18cm, integration time f = 0.1 sec, 4 
December, 1999). The baseline projection was directed mainly 
along velocity of irregularities transfering. The conditions of 
it's observations corresponded to the case of weak phase fluctu- 
ations. Spectral index p, evaluated on the slope of power spec- 
trum, is equal to 3.8. 

As stated above, proposed method theoretically may be 
used for calculation of solar wind velocity. The expression (|9j 
shows, that must be oscillations on spectrum wings. They de- 
pend on relation of velocity and baseline. Experimental spec- 
trum does not demonstrate distinct oscillations. Supposed, ve- 
locity fluctuations at solar wind were significant, and provoked 
the smoothing of minimums and rendered impossible to evalu- 
ation of velocity. The break point of measured spectrum could 
not be clear distinguished at this graph. Further theoretical and 
experimental testing is need for adjusting of procedure for solar 
wind velocity evaluation. 
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4. Conclusions 

We performed the theoretical analysis of interferometer re- 
sponse on radio emission propagated through the coronal 
plasma. The power spectrum of field fluctuations was obtained 
using the approach of "frozen-in" hypothesis for two cases: 
strong and weak phase fluctuation. It was showed that data on 
the propagation medium can be received at sounding of media 
by monochromatic signals from spacecrafts, as well as wide- 
band emission from natural radio sources. 

The theoretical analysis has demonstrated that interferom- 
eter's signal carries information about spatial spectrum of elec- 
tron density distribution, intensity of phase fluctuations and so- 
lar wind velocity. 



